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Abstract 

The dynamics of charged particles in general linear focusing lattices with quadrupole, skew- 
quadrupole, dipole, and solenoidal components, as well as torsion of the fiducial orbit and variation 
of beam energy is parameterized using a generalized Courant-Snyder (CS) theory, which extends 
the original CS theory for one degree of freedom to higher dimensions. The envelope function is 
generalized into an envelope matrix, and the phase advance is generalized into a 4D symplectic 
rotation, or an U{2) element. The ID envelope equation, also known as the Ermakov-Milne-Pinney 
equation in quantum mechanics, is generalized to an envelope matrix equation in higher dimen¬ 
sions. Other components of the original CS theory, such as the transfer matrix, Twiss functions, 
and CS invariant (also known as the Lewis invariant) all have their counterparts, with remarkably 
similar expressions, in the generalized theory. The gauge group structure of the generalized theory 
is analyzed. By fixing the gauge freedom with a desired symmetry, the generalized CS parame¬ 
terization assumes the form of the modified Iwasawa decomposition, whose importance in phase 
space optics and phase space quantum mechanics has been recently realized. This gauge fixing 
also symmetrizes the generalized envelope equation and express the theory using only the gener¬ 
alized Twiss function (3. The generalized phase advance completely determines the spectral and 
structural stability properties of a general focusing lattice. For structural stability, the generalized 
CS theory enables application of the Krein-Moser theory to greatly simplify the stability analysis. 
The generalized CS theory provides an effective tool to study coupled dynamics and to discover 
more optimized lattice design in the larger parameter space of general focusing lattices. 

PACS numbers: 29.27.-a,52.20.Dq 
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I. INTRODUCTION 


In accelerators and storage rings, charged particles are confined transversely by electro¬ 
magnetic focusing lattices. Many different kinds of focusing lattice have been successfully 
designed and implemented. The fundamental theoretical tool in designing an uncoupled 
quadrupole lattice is the Courant-Snyder (CS) theory |l|, which can be summarized as fol¬ 
lows. For a given set of focusing lattice in the x— and y—directions Kx{.t) and Hyit), particle’s 
dynamics is governed by the oscillation equation 


q + Kg{t)q = 0 , 


( 1 ) 


where q represents one of the transverse coordinates, either x or y. Solution of Eq. ([T]) can 
be expressed as a symplectic linear map M{t) that advances the phase space coordinates 




( 2 ) 


In CS theory, the linear map M{t) is given as 
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where a (t) and (3 (t) are two of the so-called Twiss parameters, and 0 (t) is the phase 
advance. They are dehned by an envelope function w (t) as 


13 (t) = (t) , 

a (t) = —ww , 

, , , dt 

Cj){t) = 


Jo P (t) 

and the envelope function w {t) is determined by the envelope equation 


(4) 

(5) 

( 6 ) 


lil + Kg (t) W = W 


-3 


(7) 


In Eq. (^, qo = q (t = 0), (jo = q (t = 0), /3o = (3 (t = 0), and ao = a (t = 0) are the initial 
conditions at t = 0. 
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Associated with the dynamics of Eq. ([T]), there exists a constant of motion, Ics, known 
as the Conrant-Synder invariant 
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Here 7 {t) is the third Twiss parameter. It turns out that the transfer matrix M (t) can be 
decomposed into the elegant form 


( 10 ) 


which seems to indicate a certain structure for M{t). 

The CS theory can be viewed as a parameterization method of the time-dependent 2 x 2 
symplectic matrix M{t) for a standard uncoupled lattice. Not surprisingly, there exist other 
parameterization schemes mathematically. Why is the CS parameterization preferable? This 
is because it describes the physics of charged particle dynamics. The main components of 
the CS theory, i.e., the phase advance, the envelope equation, the transfer matrix, and the 
CS invariant are physical quantities describing the dynamics of the particles. For example, 
the CS invariant dehnes the emittance in phase space, and the envelope function describes 
the transverse dimensions in conhguration space. This theoretical framework also makes it 
possible to investigate collective effects associated with hi^-intensity beams, such as in the 
construction of the Kapchinskij-Vladimirskij distribution [sHsJ. 

However, the CS theory can only be applied to the x— or ?/—dynamics separately for the 
ideal case of uncoupled quadrupole focusing lattices. In realistic accelerators, there exist 
bending magnets, torsion of the design orbit (hducial orbit), and skew-quadrupole compo¬ 
nents, which are introduced intentionally or by misalignment g, l7|. Solenoidal magnets 
are also used in certain applications Sj. When these additional components are included, 
the transverse dynamics in the x— and ?/—directions are coupled, and the focusing force 
depends on the transverse momentum as well. In this most general case, the transfer matrix 
M{t) is a time-dependent 4x4 symplectic matrix, which has 10 time-dependent parameters 
and admits many different schemes for parameterization. The hrst set of parameterization 
schemes for M (t) were developed by Teng and Edwards Kil and Ripken 
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which have been adopted in lattice design and particle tracking codes, such as the MAD code 
3, 3 • ^ class of different parameterizations by directly generali zing the Twiss parameters 


to higher dimensions has also been developed by Dattoli, et al. 17h19I|. However, in contrast 


to the original CS theory, these parameterization schemes are designed from mathematical 
considerations, and fail to connect with physical parameters of the beam. The elegant and 
much-needed connection with the physics of beam dynamics in the original CS theory for 
one degree of freedom is not transparent in these parameterization schemes. This is probably 
why there is no de facto standard yet adopted by the accelerator community. Another main 
reason is that for most present-day accelerators and rings, the transverse dynamics are so 
nearly decoupled that perturbative treatment often works satisfactorily. Even for lattices 
with strong coupling, elementary methods can be used to analyze the dynamics, even though 
the calculation often becomes rather involved and requires diligence and patience. 

In a recent Letter [^, we reported the development of a generalized CS theory for 
focusing lattices with the most general form in Eq. (ITTll . including bending magnets, torsion 
of the design orbit, and solenoidal magnets, in addition to quadrupole and skew-quadrupole 
components. In this generalized theory, the physics elements of the original CS theory, 
i.e., the phase advance, the envelope equation, the transfer matrix, and the CS invariant 
are all generalized to the 2D coupled case with identical structure. This new development 
also generalizes our previous results for coupled dynamics including only a skew-quadrupole 
lattice component 2ll-l24|. In this paper, we give a detailed derivation of the generalized 
CS theory reported in Ref. 20|], describe the theoretical structure of the theory in terms of 
gauge freedoms and group decomposition, and demonstrate the application of the theory in 
stability analysis. 


II. THEORETICAL MODEL AND SUMMARY OF RESULTS 

In this section, we outline the theoretical methods used and summarize the main results 
obtained in this paper. As discussed in Sec. [H when realistic components such as skew- 
quadrupoles, bending magnets, torsion of the design orbit, solenoidal magnets are included, 
in addition to the standard quadrupole components, the transverse dynamics in the x— and 
?/—directions are coupled, and the focusing force depends on the transverse momentum. In 
this case, the linear dynamics of a charged particle relative to the hducial orbit are governed 
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by a general time-dependent Hamiltonian 2^ of the form 


H = ^ ^ 


K (t) R (t) 
R (t)'^ m~^ (t) 


( 11 ) 


Here, = {x,y,Px,Py)'^ are the transverse phase space coordinates, and K{t), R(t) and 
m~^ {t) are time-dependent 2x2 matrices. The matrices A, K{t) and m~^ [t) are also 
symmetric. In this most general Hamiltonian, the skew-quadrupole and dipole components 
are included in the off-diagonal terms of the k {t) matrix, and the solenoidal component 
and the torsion of the hducial orbit are included in the R {t) matrix. There are several 


different methods to include the effect of torsion, which were reviewed by Hoffstaetter 


26| 


Typically, Frenet-Serret coordinates along the hducial orbit are used. When the hducial 
orbit is straight, the Frenet-Serret coordinates are not uniquely dehned. In this case, we 
can choose any particular set of Frenet-Serret coordinates in the straight section, as long it 
is smoothly connected to those in the curved sections. The variation of beam energy along 
the hducial orbit is rehected in the mass matrix m~^ {t ), which is allowed to be any real 
symmetric matrix for complete generality. The transfer matrix M{t) corresponding to H is 
a time-dependent 4x4 symplectic matrix, which has 10 time-dependent parameters. Our 
goal is to develop a generalized Courant-Snyder parameterization method for M{t), which 
has the same elegant structure and direct connection to beam dynamics as the original 
Courant-Snyder theory for one degree of freedom. 

We will use a time-dependent symplectic transformation technique 0, Q, Q to analyze 


the charged particle dynamics governed by the Hamiltonian given in Eq. (1111) . This technique 
is described in Sec. m The concept of scalar envelope function is generalized to a 2 x 2 
envelope matrix, and the envelope equation in 2 x 2 matrix form is developed [see Eq. fl36p ]. 
In the original CS theory, the envelope equation ([7]) is one dimensional and plays a central 


role. It also has been discovered or re-discovered many times 28l432l| in other branches of 


physics. In quantum physics, it is known as the Ermakov-Milne-Pinney 


which has been utilized to study ID time-dependent quantum systems [33|,l3J] and associated 


equation 


28|-l30|, 


non-adiabatic Berry p, 
can be found in Ref. 


liases 


35| . A brief account of the history of the ID envelope equation 


36[ |. We expect the generalization of the envelope equation to higher 


dimensions for the most general Hamiltonian to have applications in areas other than beam 
physics as well. The ID CS invariant given by Eq. (0), also known as the Lewis invariant 
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3ll. l32| in quantum physics, is generalized to higher dimensions in Eqs. flU]) and (HS]) . 

Also in Sec. lIVl the ID phase advance is generalized to a time-dependent matrix P, which 
belongs to the symplectic rotation group S'p(4) f] 5'0(4) = U{2). Here, *S'p(4), S'0(4), and 
U{2) denote the groups of 4 x 4 symplectic matrices, 4x4 rotation matrices, and 2x2 
unitary matrices, respectively. For dynamics with one degree of freedom, the phase advance 
is naturally an angle (an element of SO{2)) in the 2D phase space. For dynamics with two 
degrees of freedom, the phase space is 4D, and it is tempting to represent the phase advance 
by two angles. This is what has been adopted in previous parameterization schemes. From 
the viewpoint of theoretical physics and geometry, however, it is more natural to represent 
the phase advance for dynamics with two degrees of freedom by a 4D rotation (an element 
of S'0(4)), which is not equivalent to two 2D rotations. Because of the symplectic nature of 
the Hamiltonian dynamics, the generalized phase advance in higher dimensions thus belongs 
to the symplectic rotation group. Of course, one can adopt different views on this. In the 
normal form analysis of accelerator rings, the 4D transfer matrix is reduced to a 2D rotation 
after block diagonalization. In a sense, we can compare coupled betatron motion to the 
Dirac equation. In a fully quantum mechanical limit, the only correct approach is to treat 
electrons and positrons as inextricably coupled; only the 4D approach is permissible in this 
limit. But in the accelerator physics, we routinely treat electrons and positrons as completely 
separate entities, and two 2D descriptions are adopted without much hesitation. 

The generalized decomposition for the symplectic map M{t) is given by Eq. (03]), which 
has exactly the same structure as the original ID CS theory given by Eq. (1T0|1 . In addition 
to its aesthetic elegance, the generalized CS theory provides an effective tool to describe 
the beam dynamics governed by the most general Hamiltonian. The 2x2 envelope matrix 
w dehnes the transverse dimension of the beam, and the generalized CS invariant dehnes 
the emittance. These components of the generalized CS theory are derived in detail in 
Sec. IIVI For the present application to beam transverse dynamics, there are two degrees of 
freedom. But the theory developed is valid for any degree of freedom. For a system with 
n-degrees of freedom, the time-dependent matrix A{t) specifying the Hamiltonian in Eq. fllip 
will be 2n x 2n, the envelope matrix will be n x n, and the phase advance will belong to 
Sp{2n)r[SO{2n) = U{n). 

In Sec. E we investigate the group structure of the generalized CS theory, which is 
built on the decomposition of the time-dependent symplectic coordinate transformation G 
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in the form of Eq. 0421) . There exists a gauge freedom in this decomposition specihed by a 
2D rotation element c G 50(2) for every t. The transfer map M{t) is independent of this 
gauge. By hxing the gauge freedom with a desired symmetry, the decomposition of G as 
PS assumes the form of the modihed Iwasawa decomposition (or pre-Iwasawa decomposi¬ 
tion), whose importance in phase space optics 37, ^ and phase space quantum mechanics 


39l | has been recently realized. This specihc gauge hxing also symmetrizes the generalized 


envelope equation and express the theory using only the generalized Twiss function jS. For 
a symplectic matrix, the modihed Iwasawa decomposition is equivalent to the well-known 


Iwasawa decomposition for a semi-simple Lie group 


40|. However, the unique feature of the 


theory described here is that the decomposition is constructed as a function of time, and 
from the viewpoint of dynamics using the generalized envelope equation. Nevertheless, it is 


a pleasant surprise to hnd the deep connection between the original CS t 


particle dynamics [l| and the Iwasawa decomposition for Lie groups [^, two theoretical 
formalisms developed concurrently. This connection also demonstrates that beam dynam¬ 
ics, phase space optics and quantum dynamics have a similar theoretical structure at the 
fundamental level. In order to satisfy the symmetry requirement of the modihed Iwasawa 
decomposition, the gauge freedom need to be selected locally as a function of time, which 
is the characteristics of gauge theories in theoretical physics. This procedure also results in 
a symmetrized envelope equation in terms of the generalized Twiss parameter (3, which is a 
symmetric, positive-dehnite matrix. The beam dimensions and emittance can be expressed 
using the (3 matrix only. 

We show in Sec. IVTIhow the generalized CS theory can be used to analyze the stability of a 
charge particle dynamics in realistic accelerators with quadrupole, skew-quadrupole, dipole, 
and solenoidal components, as well as torsion of the hducial orbit and variation of beam 
energy. It turns out that the generalized phase advance as a symplectic rotation completely 
determines the spectral and structural stability properties of the general lattice after a 
matched solution of the envelope equation is found. For structural stability, the generalized 
CS theory enables us to apply the Krein-Moser theory 4ll445l| to greatly simplify the stability 
analysis. This general result includes the well-known stability criterion for sum/difference 
resonances for uncoupled quadrupole lattices as a special case. 


Ed 


leory for charged 
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III. METHOD OF TIME-DEPENDENT CANONICAL TRANSFORMATION 


We will construct the generalized Courant-Snyder theory for the general focusing lattice 
given by Eq. fllll) using a method of time-dependent canonical coordinate transformation. 
Let’s consider a linear, time-dependent Hamiltonian system with n-degrees of freedom 

H = ^z'^A{t)z, (12) 

2 ; = {Xi,X2, ...,Xn,Pl,P2, • 


Here, A (t) is a 2n x 2n time-dependent, symmetric matrix. The Hamiltonian in Eq. f fllll) i 
has this form with n = 2. The basic idea is to introduce a time-dependent linear canonical 
transformation 271 

z = S (t) z , (13) 


such that in the new coordinates z, the transformed Hamiltonian has the desired form 

H = ^z^A {t) z , 


(14) 


where A [t) is a targeted symmetric matrix. Because the transformation 0131) is canonical, 
it requires that 

SJS^ = J, (15) 


Here, J the 2n x 2n unit symplectic matrix of order 2n, 



(16) 


and I is the n x n unit matrix. Equation 0151) implies that S' is a symplectic matrix. 
In addition, it needs to satisfy a differential equation, which can be derived as follows. 
Hamilton’s equation for x is given by 


i = JVH , 


(17) 


Using index notation, Eq. flTTl) becomes 

dH 


2 “ 1 “ J • 


(18) 
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Switching back to matrix notation, Eq. fllSp can be expressed as 


z = JAz. 


( 19 ) 


Similarly, 


I = JAz = JASz. (20) 

Meanwhile, ^ can be directly calculated from Eq. flT^ by taking a time-derivative, which 
gives 

t = Sz + Sz=(S + SJA)z. (21) 

Combining Eqs. fl2Up and fl^ip gives the differential equation for S 


S = [JAS - SJA) . 


( 22 ) 


The remarkable feature of the canonical transformation S is that it is always symplectic, 
if S is initially symplectic at t = 0. This assertion can be proved by two methods. For the 


hrst proof, we follow Leach 


27| and consider the dynamics of the matrix K = SJS'^, 


k = SJS^ + SJS^ 


\JAS - SJA) JS^ + SJ 

[-SAJ + AJS^)' 

jAsjs^ - sjs^Aj 

= 

JAK - KAJ 



(23) 


Equation fl2^ has a hxed point ai K = J. If S{t = 0) is symplectic, i.e., K {t = 0) = J, then 
K = 0 and K = J for all t, and S is symplectic for all t. A more geometric proof can be given 
from the viewpoint of the flow of S. Because A is symmetric, we have JJA — A^JJ = 0, 
which indicates that JA belongs to the Lie algebra sp (2n, R). We now show that if S is 
symplectic at a given t, then JAS belongs to the tangent space of Sp{2n,R) at S, i.e., 
JAS G TsSP (2n, R). Let’s examine the Lie group right action: S : a ^ aS for any a in 
Sp {2n, R ), and the associated tangent map 


Ts : TaSp (2n, R) ^ T^gSp (2n, R ). (24) 

It is evident that jAS is the image of the Lie algebra element jA under the tangential map 
Ts- This means that JAS is a vector tangential to the space of Sp {2n, R) at S, if S is on 
Sp (2n, R). By the same argument SJA G TsSP (2n, R) as well. Thus, the right-hand side 
of Eq. (1221) is a vector on Sp (2n, R), and the S dynamics will stay on the space of Sp{2n, R). 
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We can always choose initial conditions such that S is symplectic at t = 0, and this will 
guarantee that the time-dependent transformation specihed by Eq. 0221) is symplectic for all 

t. 


IV. GENERALIZED COURANT-SNYDER THEORY 


We now apply the technique developed in Sec. Illll to the Hamiltonian system in Eq. (ITTl) . 
Our goal is to hnd a new coordinate system where the transformed Hamiltonian vanishes. 
This idea is identical to that in Hamilton-Jacobi theory. Applications of Hamilton-Jacobi 
theory include the construction of act ion-angle variables for periodic systems ^ and Ending 


geodesic curves on an ellipsoid [^. It is often required that the variables in the Hamilton- 
Jocabi equation can be separated in order for the technique to be efiective for practical 
problems. This limits its application. Since our dynamics is linear, the new coordinate 
system can be more easily constructed using the method developed in Sec. IIIIl We will 
accomplish this goal in two steps. First, we seek a coordinate transformation z = Sz such 
that, in the z coordinates, the Hamiltonian assumes the form 


It- - I /^(^) 0 

H = -z^Az , A = ^ ^ 

2 ' 0 /i(f) 


(25) 


where /r(f) is a 2 x 2 matrix to be determined. To write Eq. in the format of 2 x 2 
blocks, we let 

'^1 52 


5 = 


53 54 


and split the diEerential equation for S, i.e., Eq. fl22p . into four matrix equations, 

51 = /i53 — SiR^ -|- S 2 H , 

5 2 = /i54 — 5im“^ -|- S 2 R , 

53 = — fJjSi — S^K , 

54 = —^182 — 53771“^ -|- Si^R. 


(26) 

(27) 

(28) 
(29) 


Including /i(f), we have five 2x2 matrices unknown. The extra freedom is introduced by 
the to-be-determined /i(f). Based on the analogy with Eq. flTUD . we choose 5*2 = 0 to remove 
the freedom. We rename 54 to be tc, i.e., w = 54 , because it will be clear later that 54 is 
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the envelope matrix. Equations 0261) - 0291) become 


S, = ^ISs-S,R^, (30) 

Si = jjiwm , (31) 

4 = -^Si - SsR^ + WK , (32) 

S's = —wm + wRm, (33) 


for matrices Si, S 3 , w and fi. Because {Si, S 2 = 0, S 3 , S^ = w) describes a curve in S'p(4), 
they are consistent with the symplectic condition SiSj — S 2 S 3 = I, i.e., Siw'^ = I, which 
implies 

Si = w~'^ . (34) 


From Eq. fl3T]) . we obtain 


/i = 


(wmw'^'j 


-1 


(35) 


It is straightforward to verify that Eq. 0301) is equivalent to another symplectic condition 
S 3 SJ = S 4 S 3 . Substituting Eqs. into Eq. fl5^ . we immediately obtain the following 

matrix differential equation for the envelope matrix w. 


A 

dt 


dw 

dt 


m — wRm 




dw 

dt 


mR^ + w (^K — RmR^^ — (w'^wmw'^'^ 


-1 


0 . 


(36) 



For 77,-degrees of freedom, the envelope matrix w will be 77 x 77 , and the generalized envelope 
equation has the same form as Eq. 0361) . 

Once w is solved for from the envelope equation, we can determine Si from Eq. 034p and 
S 3 from Eq. fl5^ . In terms of the envelope matrix w, the symplectic transformation S and 
its inverse are given by 


S 

s-^ 



{w — R) mw^ 



(37) 

(38) 
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The second step is to use another coordinate transformation z = P{t)z to transform H 
into a vanishing Hamiltonian H = 0 a,t all time, thereby rendering the dynamics trivial in 
the new coordinates. The determining equation for the transformation P{t) is 


P 


-PJA 



(39) 


As explained in Sec. IIIIl the P matrix satisfying Eq. fl59]) is symplectic because JA G sp(4). 
From /i = we know that JA is also antisymmetric, i.e., JA belongs to the Lie algebra 
so(4) of the 4D rotation group S'0(4). Thus JA G sp(4)nso(4), and P{t) is a curve in 
the group of 4D symplectic rotations, i.e., P{t) G S'p(4) f] 5'0(4) = 17(2), provided the 
initial condition of P{t) is chosen such that P(0) G S'p(4) f] 5'0(4) = U{2). We call P{t) 
the generalized phase advance, an appropriate descriptor in light of the fact that P{t) is a 


symplectic rotation. The Lie algebra element (inhnitesimal generator) —JA 


0 — 
fi 0 


is 


the phase advance rate, and it is determined by the envelope matrix through Eq. 0351) . Since 
5^(4) n*S'0(4) = f/(2), P and its inverse must have the forms 


P = 


= 



-P^ 

-'2 


Pi 


(40) 


(41) 


Combining the two symplectic coordinate transformations, we obtain the transformation 


z = G{t)z = P{t)S{t)z. (42) 

In the 7 coordinate representation, because H = 0, the dynamics is trivial, i.e., 7 = const. 
This enables us to construct the symplectic matrix specifying the map between zq and 
= M{t)zo as 

w~'^ 0 \ 

[wR — w)m w j ^ 

where subscript “0” denotes initial conditions at f = 0, and Pq is taken to be I without 
loss of generality. This expression for M(t) generalizes the decomposition of the symplectic 
map for the original ID CS theory given by Eq. flTUD . The hrst and the third matrices 


M{t) = S-^p-^PoSo = 



{w — R) mw^ w 
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in Eq. (fill]) obviously have the same construction as their counterparts in Eq. The 

phase advance, as a 4D symplectic rotation in Eq. 0431) . generalizes the 2D rotation 
matrix, which is also symplectic, in Eq. OlOp . The phase advance P is generated by its 
inhnitesimal generator jA determined by the envelope matrix through /i = 

This mechanism for phase advance in 4D phase space is identical to the original ID CS 
theory where the inhnitesimal generator of the phase advance is w~‘^ for a scalar envelope 
w. The importance of the decomposition in Eqs. 043|) and flTUl) can be appreciated from 
both physical and mathematical points of view. We explain the physical meaning of the 
decomposition here, and leave the mathematical analysis to Sec. |Vl The hrst matrix from 
the right is a matching transformation at t = 0 of the initial conditions to an equivalent 
focusing system, where the phase space dynamics can be characterized by a time-dependent 
rotation. The second matrix from the right is a transformation along the time axis in this 
equivalent focusing system, with the phase advance playing the role of a time-like evolution 
parameter. And the third matrix from the right is a back-transformation to the original 
coordinate system at t > 0 . 


m 


The coordinate transformation can also be used to construct invariants of 
A general description of linear symplectic invariants can be found in Refs, 
constant 4x4 positive-dehnite matrix the quantity 

4 = z^S^P^^PSz 


he dynamics. 


49[ |. For any 


(44) 


is a constant of motion, since ^ = PSz is a constant of motion. The subscript “4’ in 4 
used to indicate that it is an invariant associated with 4 For the special case of ^ = /, the 
phase advance P in Eq. fl44D drops out, and 

/ 

Ics = z^ Sz = z^ z, (45) 



where a, /?, and 7 are 2 x 2 matrices dehned by 

a = w'^Ss , 
f5 = vj^w , 

7 = 5'JS'3 -I- . 


(46) 

(47) 

(48) 


Here, we have used Ics to denote this special invariant because it is the invariant that 
generalizes the CS invariant l| (or Lewis invariant 32|) for one degree of freedom in 
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Eq. ([H]). The matrices, a, /?, and 7 are the generalized Twiss parameters in higher dimensions. 
It is straightforward to verify that they satisfy 

Pj = I + (49) 


which is a familiar relationship in the original CS theory between the scalar Twiss parameters 
dehned by Eqs. (jl]), ([S]) and (jn]). The symplectic condition wSj = S^w'^ has been used in 
obtaining Eq. (09]). 

It has been demonstrated that the envelope matrix w and the invariant define the 


beam dimensions and emittance 
strong space-charge potential 


: or both low intensity beams and high intensity beams with 

3 - 

Note that we have “overloaded” the symbols “M, w, a, /?, 7 , /cs” to represent the 
same physical quantities in both the original CS theory for one degree of freedom and the 
generalized CS theory in higher dimensions without causing any confusion. It is actually 
more appropriate to do so than not, because the quantities in higher dimensions recover their 
counterparts for one degree of freedom as special cases, and the correspondence between 
them is exact. 


V. GROUP STRUCTURE OF THE GENERALIZED COURANT-SNYDER THE¬ 
ORY - ROTATION GAUGE AND MODIFIED IWASAWA DECOMPOSITION 

In Sec. IIV] we noted that initial conditions for the envelope matrix w need to satisfy 
the symplectic condition; otherwise they can be arbitrary. There are freedoms in the initial 
conditions and thus the solutions for w. But the transfer matrix M is independent from these 
freedoms, which are thus gauge freedoms. A subset of the gauge freedoms has the structure 
of the orthogonal group 0{n). For a time-independent element c G 0{n), we define the gauge 
transformation c : {w, P) 1 -^ [w, P) as 


w = cw , (50) 

0 t 

P = P . (51) 

I « 

Let’s show that the transformed w and P also satisfy Eqs. and fl5^ . respectively, and M 
is gauge invariant. Multiplying Eq. fl36l) by c from the left to obtain the governing equation 
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for w, 


(J/UJ ^ \ dxV rp ~ ( 7^\ { ~'T’ ~ ~T^\—^ 

—m — wRm I + —mR + w in — RmR 1 — (w wmw 1 = 0 , (52) 

CLL j CLL 

which is the same as Eq. (I5BD with w replaced by w. According to Eq. the /x matrix 
transforms as 

jl = . 

Equation (13^ thus can be re-expressed in the same form using P and fl as 




From Eqs. flHTl) and flHKI) . the S matrix and its inverse transform as 


(53) 


S = 



s-^ = s-^ 




(54) 

(55) 


Combining Eqs. (1501) and (IHTD . fl54)) and (155]) . we conclude that M is invariant under the 
gauge transformation c : {w, P) (tc, P), i.e., M = M. 

The 0{n) gauge group introduces an equivalent class ([tc], [P]) for the decomposition of 
M using {w,P). The dimension of this equivalent class is the dimension of M as a Sp{2n) 
group. To specify S' by w and w, 2n^ numbers are needed. To specify P E U{n), additional 
numbers are needed. The symplectic condition for S, SjSi = SfS^ brings (n^ — n)/2 
constrains on w and w, and the 0 {n) gauge freedom for the equivalent class is also (n^—n)/2. 
The dimension of the decomposition is therefore 

2n^ + r ?-{—^ -^-^— )=n{2n + l), (56) 


the same as the dimension of Sp{2n). 

According to the polar decomposition theorem, any non-degenerate square matrix X can 
be uniquely factored into an orthogonal matrix O and a symmetric, positive-dehnite matrix 
Q, i.e., X = OQ. As a matter of fact, Q = y/X'^X and O = XQ~^. Using this fact, at a 
hxed time t = p we can always choose c = such that w{ti) = cw{ti) 
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is symmetric. With this gauge, the canonical coordinate transformation at t becomes 

/ 


G = P 


-,-i 


(57) 


w " 0 

{wR — w)m w 

which is in the form of a modihed Iwasawa decomposition (or pre-Iwasawa decomposition) 


whose importance in phase space optics |37|, [SSj and phase space quantum mechanics 


^1 


have been recently realized. The modihed Iwasawa decomposition is the unique decomposi¬ 
tion of a 2n X 2n symplectic matrix G in the form of 


G = P 


(58) 


Y 0 
QY F-i 

where P G Sp{2n) f] SO{2n) = U{n) and Y is symmetric. Matrix Q is also symmetric, 


which is equivalent to the condition Y^'^QY = {QY)'^Y for 



to be symplectic. 


These facts are also true if the decomposition is alternatively dehned to be 

/ 


G = 


P. 


(59) 




Y 0 
QY 

For a symplectic matrix, the modihed Iwasawa decom pos ition is equivalent to the well-known 
Iwasawa decomposition for a semi-simple Lie group j40[ |. 

Making w symmetric a,t t = ti hxes the 0{n) gauge because w{ti) is unique according 
to the polar decomposition theorem. However, such a choice only makes w symmetric at 
t = ti- As in general gauge theories, we would like to pick a gauge such that the envelope 
matrix is symmetric for all t. To accommodate this desired symmetry, we need to modify 
the governing equations, especially the envelope equation. Let 

w{t) = c~^{t)u{t ), 

u{t) = \/J{Q = ^Jw^{t)w{t ), 

= w{t)u~^{t ), 

be the time-dependent polar decomposition of w(t). Here we use u(t) to denote this special 
w{t), which is symmetric and positive-dehnite for all t. The matrix u{t) is the “symmetrized” 
w(t) and equals the square-root of the generalized jS function. We will recast the envelope 
equation f | 36p in terms of /?, as in the original Courant-Snyder theory for one degree of 


freedom 


]|. The diherence is that the procedure here has to be carried out in matrix form. 
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Rewrite Eq. as 


ib + wg -|- wh = , 

(60) 

g = {fn — Rm -|- mR^)m~^, 

(61) 

h = {k — RmR^ — Rm — Rm)m~^. 

(62) 


We symmetrize Eq. fMJ]) by taking t(;^(Eq. flBUjl i + (Eq. to obtain a second order 

ordinary differential equation for /?, 

'$ — 2vlFw + nFwg + g^nFw + /S/i + [3 = . (63) 

It is a second-order equation for (3 because w^w and w^w can be expressed in terms of (3 
and $ as follows. First note that 


lii'^w = uDu + — uD‘^u — uDii , (64) 

iv^w = uDu + uii , (65) 

D = —cc~^. ( 66 ) 

Both ii and D in Eqs. fIM]) and fl6^ can be expressed as functions of [3 and 13. For ii, from 
the dehnition of u we obtain 

iiu + uii = $ , (67) 

whose left-hand side can be viewed as a linear operator on it associated with u, 


Lu{u)=uu + uu. ( 68 ) 

The properties of the linear operator L is discussed in the Appendix. Since u = is 
symmetric and positive-dehnite, is invertible to give 

u = L-^(ff ), (69) 

where L~^ is the inverse of L dehned in Eq. flA.3l) . 

To express D in terms of f3 and /3, we exam the symplectic condition 

wSj = Ss'w '^, (70) 

S 3 = —wm -|- wRm . (71) 
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Substituting in the polar decomposition w = c gives 


Lumu{D) = umuD + Dumu = {umu — iimu) + u{Rm — mR^)u . 


( 72 ) 


Therefore, 

D = {umu — iimu) + u{Rm — mR'^)u , (73) 

where u = y/P and ii = . 

Equation (1^ is a second equation for (3. Its solutions do not uniquely determine the 
envelope matrix w, which is not surprising considering that (3 = w'^w is a “symmetric” 
version of w. However, due to the 0{n) gauge freedom, (3 contains enough information to 
determine the transfer map M. In terms of u and c~^, 


. c-i 0 , 
5=1 \Su, 

0 c-i 


5, = 


u 


-1 


0 


{uR — Du — u)m u 


(74) 

(76) 


Even though the rotation matrix c{t) here is a function of t, the transformed phase advance 
is dehned the same way as in the case of a global gauge, i.e.. 


P = P 

^ U ^ 


c-^ 0 

0 c-i 


(76) 


What is modihed is the governing equation for P„, 


P = —P 

^ U ^ U 


I 


0 

0 


D 0 
0 D 


= {umu^^ 


-1 


(77) 

(78) 


The second term on the right-hand side of Eq. (1771) is due to the dependence on t of the local 
gauge. At last, the canonical coordinate transformation between 2 ; and 7 is 


Z = GZ = PySyZ , 


(79) 


and the transfer map is 


M{PP) = S-^P-^S.. 


uO • 


(80) 
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The symmetric decomposition furnished by Eqs. (1731) . (1751) . and flSUD is equivalent to the 
decomposition described in Sec. IIVI but it has three desirable features by comparison. The 
canonical coordinate transformation in Eq. fl79p has the modified Iwasawa format for all 
t. It comprises a curve of the modihed Iwasawa decomposition, developed from a dynamical 
point of view. The gauge freedom is removed, and the dimension of the symplectic transfer 
map is directly reflected by the dimension of the decomposition. At every t, M(t) is specihed 
by two n X n symmetric matrices (3 and /?, and a U{n) matrix P^. The dimension of M{t) 
is thus (n^ + n )/2 + {n^ + n)/2 + = n{2n + 1 ). 

Before ending this section, we emphasize that the purpose of studying the gauge freedom 
is to simplify the calculation of the symplectic map M and other lattice functions and beam 
parameters. By investigating the gauge freedom in the matrix envelope equation for w, 
we have found that we can actually bypass this gauge freedom and solve for the (3 matrix 
instead, which is symmetric and does not have the gauge freedom. From Eqs. fH3]) . (148 p and 
(H9l) . the generalized Twiss parameters a and 7 can also be expressed in terms of {3 and p. 
One important advantage of using the (3 matrix is that the symmetric matrices (3 and (3 form 
a linear space, which makes the numerical algorithms of searching for matched solutions for 
(3 much more efficient than for matched solutions for w. 

VI. STABILITY ANALYSIS - SPECTRAL STABILITY AND STRUCTURAL 
STABILITY 

The classical analysis by Courant and Snyder |l| on the instability induced by sum reso¬ 
nance for uncoupled transverse dynamics may give a wrong impression that coupling effects 
are always deleterious. The coupled dynamics can be stable or unstable depending on the 
specihc conhguration of the lattice, but certainly not more unstable than the uncoupled 



MeV beam energy. Our understanding of the stability properties of coupled dynamics has 
been limited by the theoretical tools available. In this section, we demonstrate how the 
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generalized Courant-Snyder theory can be applied to study the stability of the most general 
focusing lattice given by Eq. fllll) with weak and strong coupling components in realistic 
accelerators. 

For a thorough understanding, it is necessary to distinguish two types of linear stability (or 
instability). The hrst type is spectral stability, which means the linear dynamics is stable for 
all initial perturbations. The system is spectrally unstable if there exists an initial condition 
that grows without bond. In most contexts, the meaning of stability is that of spectral 
stability. The second type is the so-called structural stability (or strong stability). It mostly 
applies to systems that are spectrally stable. A spectrally stable system is structurally 
unstable if there is a spectrally unstable system inhnitesimally closed-by. Otherwise, the 
spectrally stable system is also structurally stable. The well-known result with respect to 
the stability properties of sum/difference resonances for uncoupled lattices refers to the 
structural stability under the influence of an inhnitesimal coupling component |l|. 

The spectral and structural stability of the transverse dynamics in a periodic focusing 
lattice is determined by its one-turn (or one-period) map M(T). The fact that M{T) is a 


symplectic matrix regulates the stability properties in a signihcant way 41 
the relevant results without presenting details of the proof. 


3 - 


We list here 


The spectral property is determined by the eigenvalues and their multiplicities. There 
are four possibilities: 

Cl) All eigenvalues are distinct and on the unit circle of the complex plane. 

C2) All eigenvalues are on the unit circle. There are repeated eigenvalues. But the 

geometric multiplicity for all eigenvalues is the same as the algebraic multiplicity. 

C3) All eigenvalues are on the unit circle. There are repeated eigenvalues with alge¬ 

braic multiplicity greater than the geometric multiplicity. 

C4) There exits at least one eigenvalue not on the unit circle. 

Cases C3 and C4 are spectrally unstable, and Cases Cl and C2 are spectrally stable. For 
Cases Cl and C2, we would like to know whether they are also structurally stable. I 
been shown that Case Cl is structurally stable using the symplectic nature of M(T) 

Case C2 needs to be sub-divided into two categories: 


It 

has 

41 

-45|. 
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C2.1) 


For all repeated eigenvalues, the corresponding eigenvectors have the same sig¬ 
natures. 


C2.2) 


There is at least one repeated eigenvalue whose eigenvectos have different signa¬ 
tures. 


According to the Krein-Moser theorem 4ll-l45[|. Case C2.1 is structurally stable and C2.2 
is structurally unstable. For an eigenvector ip of M{T), its signature is defined to be the 
sign of its self-product {ip, ip) = ip*iJip. The product between two eigenvectors ip and cp in 
general is defined to be {ip, (p) = ip*iJ(p, where '^’'denotes the complex conjugate of ip'^. 

To design a coupled lattice, it is desirable to be in Case Cl, which is both spectrally and 
structurally stable. As mentioned previously, for the general Hamiltonian given by Eq. fllll) . 
the parameter space satisfying this condition is large enough for most applications. Given 
a periodic lattice, we can search for a matched solution for /?, as in the original Courant- 
Snyder theory for one degree of freedom |l|. After a matched (3 is found, the one-turn map 
is 


M(T) = VP(T)-iAo, 


which implies that M{T) is similar to P{T)~^. Their eigenvalues and multiplicity are iden¬ 
tical. Because P{T) is a symplectic rotation, all of its eigenvalues are on the unit circle, 
automatically ruling out the unstable situation in C4. 

The phase advance P{T) also determines the structural stability of the system. To 
prove this assertion, let ip and (p are the eigenvectors of M{T). Then Soip and S^cp are the 
eigenvectors of P(T)~^, and 


{So'ip.Siitf, {i’A) > 


( 81 ) 


where use had been made of the fact Sois symmetric, i.e., SqJSq = J. Equation fISTD states 
that the signatures of eigenvectors of P{T)~^ and thus its structural stability are identical 
to that of M(T). 

These analyses lead to the important conclusion that the phase advance matrix P{T) 
completely determines both the spectral and structural stability of the general focusing 
lattices. This fact can significantly simplify the stability analysis in lattice design. For 
example, if the system is in Case C2, we only need to look at the signatures of the eigenvectors 
of P{T)~^ to know if it is structurally stable. According the Krein-Moser theorem, if the 


21 





eigenvectors for all repeated eigenvalues of P{T)~^ have the same signatures, then the system 
is structurally stable. Otherwise, it is structurally unstable. Let’s show that this conclusion 
recovers the classical results on the stability properties of sum/difference resonances for 


uncoupled quadrupo 
calculated to be 121 


e lattices as special cases. In this case, the phase advance matrix is 
2 ^ 


P(T)-^ 


cos 03 , 

0 

- sin 03 , 


0 sin 03 , 0 ^ 

cos (f)y 0 sin 0 y 
0 cos 03 , 0 


0 — sin (py 0 cos 0 j/ y 

where <pa: and 0^ are the one-turn phase advance in the x— and ?/—directions. Its four sets 
of eigenvalues, eigenvectors, and signatures are 


A3;+ = cos 03,i sin 03,, 03,+= (1, 0, i, 0)"^, 0-3;+= -1 , (82) 

A3;- = cos 03, - i sin 03,, 03,- = (1, 0, -0 0)^ , 0-3,- = 1, (83) 

Xy+= cos (j)y + i sin (t)y, 0J/+= (0,1,0,0'^, cry+ = -l, (84) 

Ay- = cos 0j^ — i sin 0y , 0y-= (0,1, 0, —i)"^ , Uy-= 1. (85) 


Resonance occurs when two or more eigenvalues collide, which has four possibilities: 

Rl) Self-resonance in the a:— direction. 03 , = nvr and A3;+ = A3;- = ±1. 

R2) Self-resonance in the ?/—direction. 0 y = nvr and Ay+ = Ay- = ±1. 

R3) Sum resonance. 03 , + 0 y = utt, A3;+ = Ay- and A3;- = Ay+. 

R4) Difference resonance. 03 , — 0 y = nvr, A3;+ = Ay+, and A3;- = Ay-. 

Case Rl is structurally unstable because ax+ and a^- are different. Case R2 is structurally 
unstable for the same reason. 

For the sum resonance at the repeated eigenvalue A3;+ = Ay-, the signatures <73;+ and cxy- 
of the corresponding eigenvectors have different signs. The sum resonance is thus structurally 
unstable. For the difference resonance at the hrst repeated eigenvalue A3;+ = Ay+, the 
corresponding eigenvectors 03,+ and 0 y+ have the same signature. This is also true at the 
second repeated eigenvalue A 3 ;- = Ay-. The difference resonance is thus structurally stable. 
These results are well-known previously l|, but recovered here as a special case of a more 
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general criterion based on the generalized phase advance and the Krein-Moser theorem 


li 


45[ |. We expect that the more general stability criterion expressed in terms of the phase 


advance matrix P{T) to be a powerful tool for future lattice design with strong coupling. 


VII. CONCLUSIONS AND FUTURE WORK 


We have presented in this paper a detailed derivation of the generalized Courant-Snyder 
theory for the most general linear focusing lattices with quadrupole, skew-quadrupole, dipole, 
and solenoidal components, as well as torsion of the hducial orbit and variation of beam 
energy. The theoretical structure of the theory in terms of gauge freedoms and group 
decomposition were described. We have also demonstrated the application of the theory 
in stability analysis for strongly and weakly coupled lattices. In addition to being more 
realistic, the most general Hamiltonian in Eq. (1111) enables a much larger parameter space 
for designing strongly coupled lattices that are spectrally and structurally stable. The 
generalized Courant-Snyder parameterization scheme developed here provides an effective 
tool to study the coupled dynamics and to discover more optimized lattice design in the 
larger parameter space of general focusing lattices. The formalism also sets the theoretical 
foundation for investigating collective phenomena in high-intensity beams, such as the self- 
consistent solutions of the Vlasov-Maxwell equations in p hase space including strong self-held 
effects that can couple the transverse dynamics [59|-l63j. 

As mentioned in Sec. lYl the theoretical framework developed is valid for linear system 
with any degree of freedom. In particular, we can apply it to the 3D coupled dynamics, which 
includes the sychrotron oscillation in RF cavities, and the linear coupling between transverse 


and longitudinal dynamics as in the recent investigations of emittance change 


M- 66 ]. In 


this case, n = 3 and the focusing matrix K{t) in Eq. fllip becomes a 3 x 3 matrix which 
describes both sychrotron and betatron oscillations as well as possible coupling between 
them. The envelop matrix tc is 3 x 3 and satishes Eq. fl5^ with R, m and k being 3x3 
matrices. The Twiss parameters a, /9, and 7 are 3x3 matrices, the symplectic matrix M is 
6 x 6 , and all the equations they satisfy are the same as in the case of 2 -degrees of freedom. 
Studies in these directions will be reported in future publications. 
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Appendix: Lx{y) 


In this Appendix, we derive the mathematical properties of the linear transformation used 
in Eqs. fl68p and fl72]) . Let A and X denote nxn matrices. For a symmetric, positive-dehnite 
matrix A, define the linear function associated with A on A as 

La{X) = AA + XA. 

We prove that La is invertible. It is enough to show that La is injective, i.e., La{X) = 0 
only if A = 0. Let A is in the kernel of La, i.e., 

La{X) = 0. (A.l) 

Since A is symmetric, the eigenvectors of A form a basis for vectors in K^. Expressed in this 
basis, Eq. flA.ip is 

XuV^Xu + XyV^Xu = 0, (A.2) 

where u and v are any pairs of vectors in the basis, and and A^, are the corresponding 
eigenvalues, respectively. Because the eigenvalues of a positive-dehnite matrix are positive, 
Eq. flA.2D is possible only when v'^Xu = 0. This proves that La is invertible. In terms of its 
components in this basis, L^^ is given as 

Ay 

where Y = La{X). 
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